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ABSTRACT 

We present direct constraints on how the formation of low-mass X-ray binary (LMXB) populations in galactic 
fields depends on stellar age. In this pilot study, we utilize Chandra and Hubble Space Telescope (HST) data 
to detect and characterize the X-ray point source populations of three nearby early-type galaxies: NGC 3115, 

3379, and 3384. The luminosity-weighted stellar ages of our sample span «3-10 Gyr. X-ray binary population 
synthesis models predict that the field LMXBs associated with younger stellar populations should be more 
numerous and luminous per unit stellar mass than older populations due to the evolution of LMXB donor 
star masses. Crucially, the combination of deep Chandra and HST observations allows us to test directly this 
prediction by identifying and removing counterparts to X-ray point sources that are unrelated to the field LMXB 
populations, including LMXBs that are formed dynamically in globular clusters. Galactic stars, and background 
AGN/galaxies. We find that the “young” early-type galaxy NGC 3384 («2-5 Gyr) has an excess of luminous 
field LMXBs (Lx Al (5-10) xlO 37 erg s _1 ) per unit /Gbaud luminosity (L K \ a proxy for stellar mass) than 
the “old” early-type galaxies NGC 3115 and 3379 (~8-l 0 Gyr), which results in a factor of ss2-3 excess of 
Lx/Lk for NGC 3384. This result is consistent with the X-ray binary population synthesis model predictions; 
however, our small galaxy sample size does not allow us to draw definitive conclusions on the evolution field 
LMXBs in general. We discuss how future surveys of larger galaxy samples that combine deep Chandra and 
HST data could provide a powerful new benchmark for calibrating X-ray binary population synthesis models. 
Subject headings: X-rays: binaries — X-rays: galaxies — galaxies: evolution — galaxies: elliptical and 
lenticular, cD — globular clusters: general - galaxies: stellar content 


1. INTRODUCTION 

The Chandra X-ray Obseri’atory (Chandra) has contributed 
greatly to our understanding of compact objects (e.g., neutron 
stars and black holes) through the study of X-ray binary popu- 
lations in normal galaxies. Through correlative studies involv- 
ing X-ray binary populations and galaxy physical properties 
(e.g., stellar mass and stellar ages), we have gained insight 
into the key factors that influence the formation and evolution 
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of compact objects, the stars involved in the accreting binary 
phase, and their associated remnants (e.g., millisecond pulsars 
and compact-object binaries). However, a comprehensive un- 
derstanding of the formation and evolution of compact objects 
will require empirical constraints on X-ray binary populations 
in galaxies at a variety of evolutionary stages. 

In the Milky Way, it has been observed that low mass 
X-ray binaries (LMXBs) are the dominant luminous accret- 
ing binary population, and are collectively a factor of wlO 
times more luminous than the high mass X-ray binaries (e.g., 
Grimm etal. 2002). With the advent of sub-arcsecond imag- 
ing with Chandra, it has now been shown that LMXB pop- 
ulations are prevalent in other nearby galaxies (e.g., Fab- 
biano 2006 and references therein). Chandra studies of X-ray 
binaries in local galaxies have found that the galaxy-wide in- 
tegrated X-ray power from LMXBs is strongly correlated with 
galaxy stellar mass (M*; e.g., O’Sullivan etal. 2001; Colbert 
etal. 2004; Gilfanov 2004; Kim & Fabbiano 2010; Lehmer 
etal. 2010; Boroson etal. 2011; Zhang etal. 2012). How- 
ever, in early-type galaxies, the RMS scatter in the Lx-M* 
relation is found to be «40%, which is a factor of a) 3-5 
times larger than typical measurement error and a factor of 
Ri2 times larger than that expected from statistical variations 
due to small number statistics (see Gilfanov 2004). This in- 
dicates that additional factors (e.g., stellar age and globular 
cluster LMXB systems) aside from M* alone influence the 
prevalence of LMXBs in galaxies. 

Population synthesis models predict that Lx per M* for 
LMXBs in galactic fields should decrease with increasing 
stellar age. This is primarily due to a decline in the typical 
masses of LMXB donor stars, and X-ray luminosities, with 
increasing stellar age (e.g., Belczynski etal. 2002, 2008; Fra- 
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Table 1 

Properties of Early-Type Galaxy Sample 



AT-band Position 












o'nooo 

<5.12000 

D 

a 

b 

PA 

^exclude 

log L' k 



Age 

SFR 

Galaxy 

(hr) 

(deg) 

(Mpc) 

(aremin) 

(aremin) 

(deg) 

(arcsec) 

(Lk,q) 

f K , exclude 

Sn 

(Gyr) 

(Mq yr 1 ) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(B) 

(9) 

(10) 

(11) 

(12) 

(13) 

NGC3115 

10 05 14.0 

-07 43 07 

9.68i0.41 

2.7 

1.1 

315 

20 

10.95 i0.04 

0.46 

2.3i0.5 

8-5!i:S“ 

<0.004 

NGC3379 

10 47 49.6 

+12 34 53 

10.57i0.55 

3.2 

2.7 

338 

10 

10.89i0.05 

0.25 

1.3i0.7 

8- 6 -L0 n 

<0.004 

NGC3384 

10 48 16.9 

+12 37 45 

11.59i0.77 

3.5 

1.8 

321 

20 

10.78i0.06 

0.47 

l.li().5 

3-2-lo* 

<0.009 


NOTE. — Col.(l): Target name. Col. (2) and (3): Right ascention and declination, respectively, based on A'-band data from the 2MASS large galaxy atlas (Jarrett 
etal. 2003). Col. (4): Distance and la error as determined by Tonry etal. (2001). Col.(5)-(7): A'-band galaxy major axis, minor axis, and position angle, 
respectively. Col. (8): Radius of a small circular region, centered on the galactic nucleus, that was excluded from further analyses due to both the presence of a 
bright optical background from the galactic stellar field and significant X-ray source confusion. Col. (9): Logarithm of the total galaxy-wide A'-band luminosity 
and la error. The la errors are dominated by the errors on the distances (see Column (4)). These luminosities have not been corrected for the excluded region 
specified in Column (8). Col. (10): Fraction of the total galaxy-wide A'-band luminosity contained within the excluded region. The final Lk value used in most 
calculations throughout this paper can be calculated as Lk = L' K X (1 -/k, exclude), where L’ K is the value provided in Column (9). For our sample, Lk and L’ K can be 
converted to stellar mass using M±/Lk ~ 0.66 Mg ,/A^ q, appropriate for all three galaxies. Col. (11): Globular cluster specific frequency 5jv from Harris (1991). 
Col. (12): Stellar age estimates as provided by “Sanchez-Blazquez etal. (2006) or ''Me Derm id etal. (2006). Col. (13): Integrated galaxy-wide star-formation rates 
inferred from Spitzer from Temi etal. (2009) or Shapiro etal. (2010). 


gos etal. 2008, 2013a, 2013b). Initial observational tests of 
these models have yielded conflicting results. In support of 
the population synthesis predictions, Kim & Fabbiano (2010) 
studied samples of seven young («2-5 Gyr) and seven old 
(«8-12 Gyr) nearby early-type galaxies using Chandra ob- 
servations that probe the Lx 10 38 erg s -1 end of the X-ray 
binary luminosity function (XLF). Their observations showed 
that the fraction of LMXBs that were more luminous than 
«5 x 10 38 erg s -1 is a factor of ss2 times higher («3<r) for 
the young early-type galaxy sample. Further support for 
the population synthesis expectation came from Lehmer et al. 
(2007), who utilized X-ray stacking analyses in the Extended 
Chandra Deep Field-South survey to find that optically faint 
early-type galaxies (Lp « 10 9 5 -10 10 Lb,q), a population ex- 
pected to have X-ray emission dominated by LMXBs, have 
mean LMXB emission that is a factor of «3-8 times larger 
at z ~ 0.4 (rs 4 Gyr ago) than today. More recently, how- 
ever, Zhang etal. (2012) constructed detailed LMXB XLFs 
for a sample of 20 early-type galaxies spanning luminosity- 
weighted stellar ages of 1.6-17 Gyr and found an increase 
in the A'-band luminosity normalized LMXB XLF normaliza- 
tion with increasing stellar age. The same conclusion can be 
drawn from the total LMXB Lx/Lk versus stellar age dia- 
grams constructed from the sample of 30 early-type galaxies 
studied by Boroson et al. (201 1) (see §6 and Figure 9a). These 
studies seemingly contradict previous results and the popula- 
tion synthesis expectation. 

Although the above studies and theoretical expectations 
seem to show apparent conflict, many of the differences be- 
tween studies can be explained by inhomogeneities in the 
galaxy samples and LMXB populations being studied. Al- 
though there has been substantial efforts on dynamic model- 
ing of globular clusters (GCs), published population synthesis 
modeling exclusively follows the evolution of LMXBs pro- 
duced in primordial binaries formed in galactic fields (here- 
after, field LMXBs) and does not make simultaneous pre- 
dictions for LMXBs formed in GCs. The GC LMXBs form 
efficiently through dynamical interactions of multiple bodies 
(e.g., Fabian etal. 1975; Hills 1976; Verbunt 1987), and their 
formation is expected to be tied to the dynamical and chem- 
ical properties of the host GC with little relation to the host 
galaxy stellar population age. In GCs, LMXBs are expected 


to form continuously as stars interact dynamically. While in 
galactic fields, LMXB formation and emission is expected to 
be directly connected to stellar evolution timescales of the 
secondary stars, lagging Gyrs behind previous star-formation 
episodes. 

Detailed studies of LMXB populations of local early- 
type galaxies have been substantially biased towards mas- 
sive galaxies, which tend to have rich GC systems (e.g., 
Bekki et al. 2006) and LMXB populations dominated by GC 
LMXBs and not field LMXBs (e.g., Irwin 2005). Going for- 
ward, it will be critical to test population synthesis models by 
studying specifically the field LMXB populations in a sample 
of passive galaxies with stellar populations that cover a large 
range of stellar ages. The combination of deep Chandra and 
Hubble Space Telescope ( HST ) imaging offers a means for 
identifying and separating LMXBs that are formed in galactic 
fields and GCs (e.g., Voss et al. 2009). At present, such a study 
has yet to be performed, primarily because the necessary deep 
Chandra and HST observations of early-type galaxies with 
relatively young to intermediate stellar ages («3-5 Gyr) had 
not yet been undertaken. 

To remedy the above limitation, we have conducted a pi- 
lot program consisting of deep Chandra and HST observa- 
tions of the nearby early-type galaxy NGC 3384, which has a 
luminosity-weighted stellar age of «3.2± 1.5 Gyr (e.g., Mc- 
Dermid etal. 2006; Sanchez-Blazquez etal. 2007). In this 
paper, we utilize our new Chandra and HST data to con- 
struct the field LMXB XLF for the “young” stellar popula- 
tion in NGC 3384. We compare the field LMXB XLF of 
NGC 3384 with those derived using deep archival Chandra 
and HST data of two “old” early-type galaxies, NGC 3115 
and NGC 3379, which have luminosity-weighted stellar ages 
of 8.5 ± 1.0 and 8. 6 ± 1.0 Gyr (Sanchez-Blazquez et al. 2006), 
respectively. The stellar ages for the whole sample are based 
on results from simple stellar population model fitting to both 
absorption lines and spectral energy distribution shapes. In 
this process, absorption-line strengths were all calibrated to 
the Lick/IDS system (e.g., Trager etal. 1998) and stellar ages 
were determined by interpolating between multidimensional 
model grids of absorption line indices and ages. All three 
galaxies have been shown to have smooth spatial distributions 
of luminosity-weighted stellar ages consistent with no signif- 
icant age variations across the galactic extents (Norris etal. 
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2006; Sanchez-Blazquez etal. 2007). 

Throughout this paper, we make use of three X-ray band- 
passes: 0.5-2 keV (soft band [SB]), 2-8 keV (hard band 
[HB]), and 0.5-8 keV (full band [FB]). Values of H 0 = 
70 kms' 1 Mpc -1 , = 0.3, and = 0.7 are adopted 
throughout this paper (e.g., Spergel et al. 2003). 

2. EARLY-TYPE GALAXY SAMPLE SELECTION 

Our primary goal is to measure and compare field LMXB 
XLFs down to a relatively faint X-ray luminosity limit for 
early-type galaxies spanning a wide range of stellar ages. We 
therefore required a sample of early-type galaxies with deep 
Chandra observations that have corresponding HST imaging 
across the whole galactic extents. 

Using a base sample of 376 nearby early-type galaxies from 
Ellis & O’Sullivan (2006), we searched the Chandra archive 
(from Cycle 12) for observations that were sufficient to detect 
X-ray point sources to a highly-complete 0.5-8 keV depth of 
L x r* 10 37 erg s -1 , an important limit above which a signif- 
icant number of LMXBs can be detected (see below) and a 
variety of LMXB donors and accretors are expected to con- 
tribute (Fragos etal. 2008). The Ellis & O’Sullivan (2006) 
sample was selected from the Lyon-Meudon Extragalactic 
Data Archive (LEDA) using morphological type ( T < 1.5; 
E-S0 Hubble types), distance (V < 9000 km s -1 ), and appar- 
ent magnitude (Bj < 13.5). To probe effectively the stellar- 
age dependent LMXB activity in early-type galaxies, we iso- 
lated galaxies that (1) contain passive stellar populations with 
little ongoing star-formation as measured by detailed Spitzer 
studies (e.g., Temi etal. 2009; Shapiro etal. 2010), (2) have 
small values of GC specific frequency, .S'v L 3, (where .S',v = 
■'Vgc 10°' 4(m ' ,+15:) , Nqc is the number of GCs within the galaxy, 
and My is the absolute V-band magnitude of the galaxy), to 
study optimally both field and GC LMXBs, but without GC 
LMXBs dominating (e.g., Kim & Fabbiano 2004; Irwin 2005; 
Juett etal. 2005; Boroson etal. 2011; Zhang etal. 2012), (3) 
have distances D 4S 15 Mpc, as to avoid serious confusion is- 
sues, and (4) have V-band luminosities (a proxy for stellar 
mass) in a relatively small range « 2-10 x 10 111 
The final V-band luminosity criterion removes large varia- 
tions in LMXB populations that might arise as a result of the 
Lx-M* relation, ensures that a sufficient, statistically mean- 
ingful number of LMXBs ( a) 50-100) per galaxy are detected 
above Lx ~ 10 37 erg s -1 , and minimizes contributions from 
hot X-ray gas that are prevalent in more massive early-type 
galaxies (David et al. 2006). 

From the criteria above, we found that three galaxies, 
Cen A, NGC 3115, and NGC 3379 had sufficient Chandra 
coverage to study in detail LMXB XLFs; however, Cen A 
is problematic for our program, since it has highly uncer- 
tain global stellar population ages due to large dust-lanes, 
shows signatures of recent gas-rich merger activity (e.g., 
Kainulainen et al. 2009), and further contains bright extended 
X-ray structures (shells and a jet) making it difficult to study 
the XRB populations (e.g., Voss etal. 2009). By contrast, 
NGC 3115 and 3379 have deep Chandra and HST expo- 
sures and satisfy all four selection criteria posed above and 
are therefore ideal candidates to include in our program (see, 
e.g., Kim etal. 2006, 2009; Brassington etal. 2008, 2010). 

15 AXband luminosities were computed using data from the Two Micron 
All Sky Survey (2MASS) large galaxy atlas (Jarrett etal. 2003) and AXband 
extended source catalogs. The 2MASS extended source catalog is available 
online via http://tdc-www.harvard.edu/catalogs/tmxsc.html 


These galaxies have been determined to have relatively old 
stellar populations («8-10 Gyr based on multiple stellar age 
measurements) and therefore these Chandra and HST obser- 
vations have the potential to provide key observational con- 
straints for “old” field LMXB populations. 

Using the remaining early-type galaxy sample that did not 
have deep Chandra data in the archive, we searched for 
“young” ( ,$ 5 Gyr) early-type galaxy candidates to target 
with Chandra and HST. We found that NGC 3377 and 3384 
were the nearest sources (D < 15 Mpc) that have relatively 
young stellar populations and satisfy our selection criteria. 
In Cycle 13, we successfully proposed for joint Chandra and 
HST observations of NGC 3384, which we present for the first 
time in this paper. Taken together, NGC 31 15, 3379, and 3384 
form a clean sample of passive D < 15 Mpc early-type galax- 
ies with stellar ages ranging from rG- 10 Gyr. The Chandra 
and HST observations of all three galaxies cover the entire V- 
band extents and reach 0.5-8 keV limits of 10 37 erg s -1 . 
The V-band luminosities provide good proxies for the stellar 
masses M* of the galaxies in our study. Adopting the prescrip- 
tion outlined in Appendix 2 of Bell et al. (2003), we utilized 
B — V colors to estimate that M k / L K « 0.66 Mq/L^© for all 
three galaxies in our sample. We note that this method has 
an uncertainty of «0. 1-0.2 dex on the derived M+/Lk for a 
given B — V color. It is likely the case that the younger galaxy 
NGC 3384 has a true M+/Lk that is lower than those for the 
two older galaxies NGC 3115 and 3379. Hereafter, we make 
use of scalings with Lk , and where appropriate, convert to M* 
using a constant M k ! L K conversion factor for the whole sam- 
ple. In Table 1, we summarize the properties of the galaxies 
that make up our sample and in § 3 we describe our Chandra 
and HST data and analyses. 

3. DATA REDUCTION 

As discussed in § 2, we are comparing results from new 
Chandra and HST observations of the “young” early-type 
galaxy NGC 3384 with comparable-to-better quality archival 
data for “old” early-type galaxies NGC 3115 and NGC 3379. 
We note that there are published papers describing some of the 
Chandra and / 1ST archival data available for NGC 3115 and 
NGC 3379 (e.g., Brassington etal. 2008, 2010, 2012; Wong 
etal. 2013); however, to limit comparison errors, we analyze 
the Chandra and HST data for all three galaxies using the 
methodologies adopted in this paper. In the sections below, 
we describe our data reduction methods in detail. 

3.1. Chandra Data Reduction and X-ray Catalog Creation 

The Chandra observations for all three galaxies were con- 
ducted in multiple observations (hereafter, ObsIDs) using the 
S-array of the Advanced CCD Imaging Spectrometer (ACIS- 
S; see Table 2 for full observation log). The V-band semi- 
major axes of the galaxies (Column 5 of Table 1) range from 
2.7 , -3.5' and therefore subtend an angular extent smaller 
than a single ACIS-S chip («8.'5 x 8(5). As such, for all 
Chandra observations, the galactic extents were almost en- 
tirely contained within the ACIS-S3 chip for all observations. 
For our data reduction, we made use of Cl AO v. 4.5 with 
CALDB v. 4.5.6. We began by reprocessing the pipeline pro- 
duced events lists, bringing level 1 to level 2 using the script 
chandra_repro. The chandra_repro script runs a va- 
riety of Cl AO tools that identify and remove events from 
bad pixels and columns, and filter the events list to include 
only good time intervals without significant flares and non- 
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Table 2 

Chandra Advanced CCD Imaging Spectrometer (ACS) Observation Log 


Galaxy 

Obs. ID 

Obs. Start 
(UT) 

Exposure Time" 
(ks) 

Flaring 6 

Intervals 

Flaring Time 6 
(ks) 

Obs. Mode c 

NGC3115 

2040 

2001 Jun 14, 10:31 

35.3 

1 

1.6 

F 


11268 

2010 Jan 27, 00:42 

40.6 

0 


VF 


12095 

2010 Jan 29, 14:29 

75.7 

0 


VF 


13817 

2012 Jan 18, 14:08 

172.0 

0 


VF 


13819 

2012 Jan 26, 22:21 

70.0 

2 

5.5 

VF 


13820“' 

2012 Jan 31, 14:55 

184.2 

0 


VF 


13821 

2012 Feb 3,09:39 

157.5 

1 

0.5 

VF 


13822 

2012 Jan 21, 08:54 

143.0 

2 

17.2 

VF 


14383 

2012 Apr 4,02:37 

119.0 

1 

0.5 

VF 


14384 

2012 Apr 6,17:51 

69.7 

0 


VF 


14419 

2012 Apr 5,18:49 

46.3 

0 


VF 


Merged 


1113.3 

7 

25.3 


NGC 3379 

1587 

2001 Feb 13,01:33 

30.5 

2 

1.0 

F 


7073 d 

2006 Jan 1,05:50 

81.6 

3 

2.5 

F 


7074 

2006 Apr 9, 10:34 

69.1 

0 


VF 


7075 

2006 Jul 3,11:04 

82.1 

2 

1.0 

VF 


7076 

2007 Jan 10,01:51 

69.2 

0 


VF 


Merged 


333.5 

7 

4.5 


NGC 3384 

4692 

2004 Oct 19, 18:57 

9.9 

0 


VF 


11782 

2010 Jan 19, 11:31 

28.2 

1 

0.5 

F 


13829“' 

2012 Jun 23, 11:07 

94.8 

0 


VF 


Merged 


132.9 

1 

0.5 



Note. — Links to the data sets in this table have been provided in the electronic edition. 

a All observations were continuous. These times have been corrected for removed data that was affected by high background; see § 3.1. 
b Number of flaring intervals and their combined duration. These intervals were rejected from further analyses. 
c The observing mode (F=Faint mode; VF=Very Faint mode). 

d Indicates Obs. ID by which all other observations are reprojected to for alignment purposes. This Obs. ID was chosen for reprojection as it had the longest 
initial exposure time, before flaring intervals were removed. 


Table 3 

HST Advanced Camera for Surveys (ACS) Observation Log 


Line Dither Pattern 


Galaxy 

Field 

Number 

Obs. Start 
(UT) 

Exposure Time (s) 
(F475W) (F850LP) 

(F475W) 

V pts 

(F850LP) 

Spacing 

(arcsec) 

NGC 31 15 

1 

2012 Mar 7, 02:49 

824 

1170 

2 

3 

3.011 


2 

2012 Mar 7,04:25 

824 

1170 

2 

3 

3.011 


3 

2012 Mar 7, 06:14 

722 

1137 

2 

4 

3.011 


4 

2012 Mar 7, 08:06 

824 

1170 

2 

3 

3.011 


5 

2012 Mar 9, 17:05 

824 

1170 

2 

3 

3.011 


6 

2012 Mar 10, 00:07 

824 

1170 

2 

3 

3.011 

NGC 3379 

2 

2010 Mar 25, 04:08 

760 

1285 

2 

3 

0.145 


3 

2010 Apr 01, 00:38 

754 

1275 

2 

3 

0.145 


4 

2010 Apr 01, 02:14 

754 

1275 

2 

3 

0.145 


5 

2010 Apr 01, 03:50 

754 

1275 

2 

3 

0.145 


6 

2010 Apr 01, 05:26 

754 

1275 

2 

3 

0.145 

NGC 3384 

4 

2012 Apr 10, 14:30 

780 

1305 

2 

3 

2.994 


5 

2012 Apr 18, 08:59 

800 

1290 

2 

3 

2.994 


6 

2012 Apr 14, 03:06 

780 

1305 

2 

3 

2.994 


Note. — Links to the data sets in this table have been provided in the electronic edition. 
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NGC 3115 (1.1 Ms) 

• 

• • • •% 


2-4 keV 

•• " . 

4-*8 keV 




NGC ^384 (130 ks) 


• * . 


. V 


\ 

• i 

2-4 keV N * 

4-8 keV 



Figure 1. Three-color Chandra images of the three early-type galaxies in our sample. Each image was constructed from 0.5-2 keV {red), 2-A keV {green), 
and 4-8 keV {blue) exposure-corrected adaptively smoothed images. Dashed ellipses represent the total A^-band galaxy size and orientation as described in 
Columns (2)-(3) and (5)-(7) of Table 1. We note that the range of the Chandra image depths (see annotations) is large (0. 1-1.1 Ms), thus affecting comparitive 
quality. 




Figure 2. Color HST ACS mosaics of the three early-type galaxies in our sample. Each image was constructed using F850LP {red), (F850LP + F475W)/2 
{green), and F475W {blue) image mosaics. Image sizes, orientations, and dashed ellipses are the same as in Figure 1. Each 202" x 202" ACS field of view has 
been outlined {solid squares) and the field name annotated. A description of our mosaicking procedure can be found in § 3.2. 
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cosmic ray events corresponding to the standard ASCA grade 
set ( ASCA grades 0, 2, 3, 4, 6). 

Using the reprocessed level 2 events lists for each Ob- 
sID, we generated preliminary FB images and point-spread 
function (PSF) maps (using the tool mkpsfmap) with a 
monochromatic energy of 1 .497 keV and an encircled counts 
fraction (ECF) set to 0.393. For each ObsID, we created pre- 
liminary source catalogs by searching the FB images with 
wavdetect (run including our PSF map), which was set at 
a false-positive probability threshold of 1 x 10 -5 and run over 
seven wavelet scales from 1-8 pixels (1, y/2, 2, 2\fl, 4, 4\/2, 
and 8). To measure sensitively whether any significant flares 
remained in our observations, we constructed point-source- 
excluded 0.5-8 keV background light curves for each ObsID 
with 500 s time bins. For each galaxy, we found 1-7 inter- 
vals across all ObsIDs with flaring events of gi, 3 a above the 
nominal background; these intervals were removed from fur- 
ther analyses (see Table 2). In general, this resulted in the 
removal of short intervals; however, for NGC 3115 ~25 ks of 
data were excluded. The resulting cumulative exposure times 
for NGC 3115, 3379, and 3384 were 1,113 ks, 333 ks, and 
133 ks, respectively. 

For each galaxy and ObsID, we utilized the preliminary 
source catalogs based on wavdetect positions to register 
each good-time-interval screened aspect solution and events 
list to the frame with the longest exposure time (see Table 2); 
this was achieved using Cl AO tools repro ject_aspect 
and repro ject_events, respectively. The resulting as- 
trometric reprojections gave nearly negligible linear trans- 
lations (<0.38 pixels), rotations (<0.1 deg), and stretches 
(<0.06% of the pixel size) for all ObsIDs. We created merged 
events lists for each galaxy using the reprojected events lists 
and the Cl AO script merge_obs. 

Using the individual ObsID and merged events lists, we 
constructed images for each galaxy in the FB, SB, and HB. 
For each of the three bands, we created corresponding ex- 
posure maps following the basic procedure outlined in §3.2 
of Hornschemeier etal. (2001); these maps were normal- 
ized to the effective exposures of sources located at the aim 
points. This procedure takes into account the effects of vi- 
gnetting, gaps between the CCDs, bad column and pixel fil- 
tering, and the spatially-dependent degredation of the ACIS 
optical blocking filter. A photon index of T = 1 .4, the slope 
of the extragalactic cosmic X-ray background in the full band 
(e.g., Hickox & Markevitch 2006), was assumed in creating 
the exposure maps. 

For a given galaxy, we created a candidate source catalog 
by first searching each image (i.e., all ObsIDs and merged 
images in all three bands) using wavdetect at a liberal 
false-positive probability threshold of 1 x 10 -5 over the s/2 
sequence (see above). We searched all ObsIDs to isolate the 
locations of transient source candidates that may have been di- 
luted by background and thus not detected by wavdetect in 
the merged image. We then merged the wavdetect catalogs 
together by cross-matching the point-source catalogs. For a 
given galaxy, we cross-matched the catalogs using matching 
radii of 1 "5, 2"5, and 3"5 for sources offset by £ 2', 2'-6', 
and 6', respectively, from the exposure-weighted average 
aim-point location. We inspected all the images and candi- 
date source regions by eye to see if any additional source can- 
didates were missed. In particular, in the crowded nuclear 
regions of the elliptical galaxies, low-flux point sources may 
not be picked up by wavdetect due to the complex back- 


grounds. As a result of our inspection, we identified six addi- 
tional candidate sources, including 3, 1, and 2 sources in the 
nuclear regions of NGC 3115, 3379, and 3384, respectively, 
which we added to our candidate source catalogs. 

We performed detailed photometry on our candidate 
source catalogs using the ACIS EXTRACT (hereafter, AE; 
Broos etal. 2012) point-source analysis software and the 
wavdetect-based positions^ The ae software contains al- 
gorithms that allow for appropriate computation of source 
properties when multiple observations with different roll an- 
gles and/or aim points are being combined and analyzed (see 
discussion below for further details), ae also uses complex al- 
gorithms for accurately computing the photometry of sources 
in crowded regions like those found in the nuclear regions of 
our early-type galaxies. The improved photometry provided 
by ae allows us to evaluate clearly the significances of the point 
sources in each candidate source catalog and prune our source 
lists to include reliable sources. For a source to be considered 
detected, it must contain s counts within an aperture repre- 
senting «90% of the point-source encircled-energy fraction 
(EEF) satisfying the following binomial probability criterion: 


P (x >s) = 


z -— ' x\(n—x)\ 

X=S 


pV - P r X < thresh, 


( 1 ) 


where n = s + /? ext and P = 1/(1 +^ext/^src)- Here b cxl is the 
total number of background counts extracted from a large re- 
gion outside of the point source (while masking out regions 
from other X-ray detected sources) that was used to obtain 
an estimate of the local background count rate. The quan- 
tity b SIC is the estimated number of background counts within 
the source extraction region, which was measured by rescal- 
ing b ext to the area of the source aperture (the inter-quartile 
range of b ext /b SIC « 24-118). We adopted a threshold value 
of P thresh = 0.004, below which we considered a source can- 
didate to be a valid X-ray point source and unlikely to be 
a fluctuation of the background. Our choice was motivated 
by the Xue etal. (2011) analysis of the 4 Ms Chandra Deep 
Field-South survey, in which the use of //hresh = 0.004 max- 
imized the number of valid sources that had optical/near-IR 
counterparts without introducing a significant number of false 
sources. 

For each galaxy, we constructed main catalogs by iter- 
atively pruning our candidate source lists to include only 
sources with P < P t hresh in at least one of the three bandpasses. 
In this process, we first computed P using only the merged 
photometry^] we then re-extracted source and background 
counts from the remaining sources, and repeated the process 
until no further sources were removed. The three galaxies re- 
quired 1-3 iterations before converging. 

We note that our multi-stage procedure for identifying a 
highly reliable list of source candidates will not result in a 
complete selection of all real sources in each image having 
binomial probabilities < fVesh, since the initial wavdetect 
selection has more complex source-detection criteria (see 
Freeman et al. 2002 for details) than the simple criterion given 


16 The ACIS EXTRACT software is available on the WWW at 
http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.html 

1 ' Had we chosen to consider P computed from both the merged observa- 
tions and each of the individual observations, additional faint variable sources 
would be included in our main catalogs. However, the deliberate inclusion of 
such sources would result in biasing our XLFs upwards for increasing num- 
bers of ObsIDs. To some extent, bright variable sources may have a small 
effect on our results, if the sources are detected in the merged images. 
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in equation 1. Such a caveat is important when considering 
completeness of our X-ray catalogs. In § 5 we discuss our 
methods for correcting for these completeness issues when 
computing LMXB XLFs. 

3.2. HST Data Reduction and Optical Catalog Creation 

Each of the three galaxies was observed with HST using the 
Advanced Camera for Surveys (ACS) with the F475W and 
F850LP filters; hereafter the £475 and zsso bands, respectively. 
The ACS field of view subtends an area of ft202" x 202"; 
therefore the full mapping of each galaxy required multiple 
pointings. For a given galaxy, each ACS pointing pattern 
was carefully chosen to cover the entire optical extent of the 
galaxy and achieve higher sensitivity near the galactic center. 
High sensitivity in the central regions of the galaxies is cru- 
cial for effectively identifying GCs and unrelated background 
candidates that fall within the bright stellar emission near the 
galactic centers. Details regarding the ACS observations are 
logged in Table 3, and in Figure 2 we show mosaicked ACS 
images of the three galaxies in our sample with each of the 
202 " x 202 " observational fields highlighted. 

For each galaxy, we first constructed a mosaicked image in 
the ^475 band by aligning each of the images to a chosen ref- 
erence observation and “drizzling” (Fruchter & Hook 2002) 
the images to the common frame. For NGC 3115, we first 
aligned ACS fields 3 and 4 (see annotations in Fig. 2) to cre- 
ate a large central field that contained some overlapping areas 
common to the remaining ACS fields 1, 2, 5, and 6 ; these 
fields were subsequently aligned to the large central field. For 
NGC 3379 and NGC 3384, we aligned all observations to 
the central images (i.e., ACS fields 2 and 5 for NGC 3379 
and NGC 3384, respectively). We performed field align- 
ments using the Pyraf tools Tweakreg and Tweakback, 
which were available through the drizzlepac v.4.3. (EJ 
software package. Tweakreg and Tweakback identify 
sources that are common to each ACS image and a chosen 
reference image (here the galactic centers) and update the 
image headers once an astrometric solution is found. Given 
the small overlaps between images, we implemented only 
small shifts in right ascension and declination to align our 
images. In general, ft4-200 bright sources per image were 
used in the field alignments, and the resulting shifts and 
RMS residuals to the fits were ^ 10 ACS pixels ( 0"5) 
and ^,0.1 pixels ( f 0"005), respectively. After aligning 
the ACS fields, we constructed #475 mosaicked images us- 
ing the astrodrizzle tool within drizzlepac. The 
astrodrizzle procedure uses the aligned, flat-field cali- 
brated and charge-transfer efficiency (CTE) corrected images 
to create a distortion-corrected mosaicked image with bad 
pixels and cosmic rays removed. 

We repeated the above ACS image mosaicking procedure 
for the Z 850 band data, with the exception that all ACS obser- 
vations of a given galaxy were aligned to the full y 475 mo- 
saicked image. The resulting two-colored mosaicked images 
are displayed in Figure 2 with the AG band galactic extents 
overlaid. 

To construct source catalogs, we first searched each im- 
age mosaic using SExtractor (Bertin & Arnouts 1996). 
We chose the following SExtractor parameters: mini- 
mum detection area of 5 pixels; detection and photometry 
analysis thresholds of 3<r; background grid of 16 x 16 pix- 
els with a filtering Gaussian of 2.5 pixels full width at half- 

18 See http://drizzlepac.stsci.edu/ for details on drizzlepac. 


Table 4 

Summary of X-ray Source Classifications in Galaxy Footprints 


Galaxy Mot AW ^ ^lmxb 

0 ) ( 2 ) ( 3 ) ( 4 ) ( 5 ) ( 6 ) 

NGC 31 15 153 46 25 9 73 

NGC 3379 151 17 10 18 106 

NGC 3384 65 18 3 14 30 


NOTE. — Col.(l): Galaxy name. Col. (2) Total number of X-ray sources de- 
tected within the galaxy footprint, including the nuclear regions. Col. (3) 
Number of nuclear sources that were excluded from XLF computations (see 
§§ 3.2 and 5.2 for explanation and justification). After excluding the nuclear 
sources, Col.(4)-(6) provide the number of X-ray detected sources classified 
as GC, unrelated background, and field LMXB candidates, respectively. 

maximum (FWHM); and 32 deblending sub-thresholds. The 
SExtractor catalogs include some non-negligible contri- 
butions from background fluctuations and unrejected cosmic 
rays (e.g., in chip gap regions where dithering makes cosmic- 
ray rejection difficult). For each galaxy, we constructed a 
master optical catalog, in which we required that sources be 
detected in both the £475 and z 850 bands, using a matching ra- 
dius of 0"5. 

Within the central 10"-20" of the galactic centers, the 
stellar field backgrounds are very high, particularly in the 
Z 850 band, and our SExtractor searching does not yield 
source detections in these regions. Since we are unable to 
identify optical counterparts to X-ray sources (e.g., GCs and 
background sources) in these regions, we reject the X-ray 
sources in these regions from consideration when computing 
field LMXB XLFs. As we will describe in §5, these regions 
are also subject to significant X-ray point source confusion, 
which supports our choice to reject these regions. The sizes 
of these rejected regions have been tabulated in Column (8) 
of Table 1, and the fraction of the X-band luminosity con- 
tained within these regions, as measured by the X-band data, 
are provided in Column (10) in Table 1. Hereafter, we utilize 
X-band luminosities that have been corrected for the exclu- 
sion of these regions. 

3.3. Astrometric Alignment of Chandra and HST Frames and 
Catalog Summaries 

We further refined the absolute astrometry of our Chan- 
dra and HST data products and catalogs by registering them 
to either the Eighth Sloan Digital Sky Survey (SDSS) Data 
Release (DR 8 ; Aihara etal. 2011) frame (NGC 3379 and 
3384) or the United States Naval Observatory CCD Astro- 
graph Catalog (UCAC2; Zacharias etal. 2000) frame when 
SDSS DR 8 data were unavailable (NGC 3115). SDSS DR 8 
was astrometrically calibrated to UCAC2 frame; however, the 
SDSS DR 8 reaches a much deeper limiting /'-band magni- 
tude of rim, » 22.5 (50% complete), compared to rn m ft 16 
for UCAC2, and is therefore our preferred optical reference 
frame. 

To align our HST frames with reference optical frames (i.e., 
SDSS DR 8 or UCAC2), we first cross-matched the HST mas- 
ter optical catalogs of each galaxy with the reference catalog 
using a 1"5 matching radius. We then applied median astro- 
metric shifts in right ascension and declination (total offsets 
ranged from 0"25-0"56 for the three galaxies) to the HST im- 
ages and source catalogs to bring them into alignment with the 
reference frame. To align the Chandra catalogs and data prod- 
ucts, we matched the Chandra main catalogs of each galaxy 
to their corresponding astrometry-corrected HST master op- 
tical catalogs using a matching radius of l."5. Again, sim- 
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Figure 3. Sample ACS image postage stamps of X-ray detected GC ( top panels) and unrelated background source (.bottom panels) candidates in NGC 3379. 
Our optical source counterpart classification criteria are described in § 4. Each postage stamp is centered on the optical source position and spans a 1 " 5 X 2”5 
region. For scaling purposes, we have included a circle with \" radius centered on the location of the X-ray source. 


pie median shifts in right ascension and declination (offsets 
ranged from 0"17-0"26 for the three galaxies) were applied 
to the Chandra data products and catalogs to bring them into 
alignment with the HST and reference optical frames (by ex- 
tension). 

The final HST and Chandra image and catalog registrations 
have a 1 a error of ft0"12. Our further analyses are focused on 
sources that lie within the X-band extents of the galaxies out- 
side the crowded nuclear regions (see Table 1). Within these 
galactic footprints, our main Chandra catalogs (see §3.1) con- 
tain 107, 134, and 47 X-ray sources for NGC 3115, 3379, 
and 3384, respectively, with the faintest FB detected sources 
reaching flux levels of 5.9, 7.9, and 16 x 10 -17 erg cm -2 s -1 , 
respectively. 

The optical limits of our HST catalogs vary substantially 
across the galactic footprints due to large background gra- 
dients from the galactic stellar fields. Within the galactic 
footprints (after excluding the inner regions specified in Ta- 
ble 1), our optical catalogs contain 582, 819, and 1125 total 
sources, which imply source densities of 64.8, 30.2, and 57.9 
arcmin -2 for NGC 3115, 3379, and 3384, respectively. As we 
will point out later, the factor of ft2 times lower source den- 
sity for NGC 3379 is also evident in the background X-ray 
number counts and is may be due to cosmic variance on 
the (small) size scale of the galactic footprints (see §5.3). 
Within the galactic footprints, we find that the #475 -band lim- 
iting depths of our catalogs is deepest in the innermost re- 
gions (outside the excluded radii given in Column (8) in Ta- 
ble 1) where the overlapping HST coverage generally con- 
tains a factor of ~2 times deeper imaging than the galactic 
outskirts. Despite the deeper imaging in the central regions, 
we find that the zsso-band limiting depths are shallowest in 
the inner regions of the galaxy due to the enhanced back- 
ground from the host galaxies in these regions. Specifically, 
in the inner regions of the galaxies, we find 5a limiting depths 
of gj™ « (26.0,26.0,26.3) and z 1 ™ ft (24.5,24.5,24.7) for 
NGC 3115, 3379, and 3384, respectively. In the outer regions, 
the respective 5er limiting depths are g 4™ ft (25.6,25.6,26.0) 
and Zgfo ft (24.9,24.9,25.1). 


4. X-RAY SOURCE CLASSIFICATION 

For each galaxy, we cross-correlated our astrometry- 
corrected Chandra main catalogs with our HST ACS optical 
master catalogs to identify likely counterparts to a subset of 
the X-ray point sources. An optical source was considered 
to be a likely counterpart if its position was located within 
the 95% confidence radius r 95 of the X-ray source. The size 
of >'95 depends on both net source counts and positional off- 
set from the exposure-weighted mean aim point; these depen- 
dencies were calculated following equation (12) of Kim etal. 
(2007). In total, we found 34, 28, and 17 optical counterparts 
to the X-ray point sources that were within the galactic extents 
of NGC 3115, 3379, and 3384, respectively; these represent 
32%, 21%, and 36% of the X-ray sources, respectively. For 
each galaxy, we estimated the total number of false matches 
expected, given our matching criteria. To do this, we first esti- 
mated the optical source density p opt as a function of angular 
distance from the galactic centers. The expected number of 
false matches was then calculated as lVf a ise ~ E,- p 0 pt,i x E^Ay, 
where p opU is the density of optical sources in the ith annu- 
lus and Ajj is the area enclosed by the 95% confidence posi- 
tional error circle of the /th source residing in the ith annu- 
lus. We estimate that /V| a | se ft l.l,2.4,and 1.3 for NGC 31 15, 
3379, and 3384, respectively. These values imply that ft6% 
of the X-ray sources with optical counterparts could be false 
matches, which will not have a material effect on our results. 

We expect that optical counterparts to the X-ray sources 
will fall in one of two broad categories: LMXBs in GCs (ex- 
pected to have zsso ~ 18-24 mag, safely above the sensitiv- 
ity limits in §3.3) and unrelated background sources. Here- 
after, we define unrelated background sources to comprise all 
X-ray sources that are not intrinsically related to the galax- 
ies being studied here. This category therefore includes fore- 
ground Galactic stars, as well as background AGN and galax- 
ies. Given that the most optically luminous counterparts to 
field LMXBs are expected to be normal stars with g 475 and 
Z850 magnitudes well below our detection thresholds (e.g., lu- 
minous red giants will have zsso AL 28-3 1 mag at the dis- 
tances to our galaxies), we do not expect to detect optical 
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counterparts to individual field LMXBs. We therefore con- 
sider X-ray sources without optical counterparts to be candi- 
date field LMXBs. 

For X-ray sources with optical counterparts, we per- 
formed basic source classifications using their optical col- 
ors and light-profile information. Optical sources were clas- 
sified as likely GCs if they had (1) colors in the range of 
0-6 < ^475— Z850 < 1.6; (2) absolute magnitudes (based on the 
distances to each galaxy) in the range of -12. 5 <M Z < -6.5; 
(3) extended light profiles in either the £475 or zsso band- 
passes, characterized as having SExtractor stellarity pa- 
rameters CLASS_STAR < 0.9 or aperture magnitude differ- 
ences m{r\) — m{r 2 ) > 0.4, where apertures consist of circles 
with radii r\ = 0."25 and r^ = 075; and (4) light profiles in 
both the £475 and -;850 bands that were not too extended to be 
GCs, defined as m{r\)—m{r 2 ) < 0.9. X-ray sources with opti- 
cal counterparts that did not meet these criteria were classified 
as unrelated background objects. 

In Figure 3, we show ten random examples of X-ray sources 
in NGC 3379 classified as GC and unrelated background 
sources. We visually inspected images similar to those shown 
in Figure 3 for all optical counterparts in all three galaxies. We 
found that our adopted selection criteria yielded reasonably 
well segregated samples of GCs and unrelated background 
sources with few sources where the classifications are ob- 
viously incorrect. Furthermore, the inferred 0.5-8 keV lu- 
minosities of the GC counterparts were all less than Lx ~ 
7 x 10 38 erg s _1 , which is well within the expected range for 
GC LMXBs. We note that it is likely that a small minority of 
sources will have incorrect source classifications (e.g., faint 
sources with poor constraints on spatial properties); however, 
we do not expect these misclassifications to have any material 
effects on our conclusions. We also note that our ability to 
identify and separate GC LMXBs from field LMXBs is de- 
pendent on the completeness of our GC catalogs. Using the 
full HST catalogs, we constructed observed optical luminos- 
ity functions (i.e., number of detected sources per magnitude 
as a function of absolute magnitude M-) for sources detected 
within the galactic footprints. For all three galaxies, we found 
steadily rising number counts with decreasing M z down to 
limits of M hm « -5.0 to -5.5 before incompleteness began 
to affect source detection. As discussed above, GCs were se- 
lected to to have -12.5 < M z < -6.5, boundaries that are well 
above the highly-complete limits. Our observed GC optical 
luminosity functions are consistent with those of the Milky 
Way, which to first order can be described as a Gaussian with 
peak luminosity at M z ss -7.9 mag and <r « 1.3 mag (Harris 
1991). We therefore expect our GC selections to be highly 
complete. 

In Table 4, we summarize our optical source categorizations 
for the X-ray sources. As discussed above, we have excluded 
from consideration sources that were within small regions of 
the galactic centers due to the very bright stellar fields in these 
regions (see §3.2). The number of excluded sources has been 
provided in Column (3) of Table 4. Within the acceptable re- 
gions of the galaxy extents, we find 30-106 field LMXBs and 
3-25 GC LMXBs per galaxy. In the sections below, we de- 
scribe our procedure for generating LMXB XLF populations 
and present our resulting XLFs. 

5. X-RAY LUMINOSITY FUNCTION COMPUTATION AND FITTING 

We chose to compute XLFs in the FB (0.5-8 keV) to make 
direct comparisons with recent XLF studies of early-type 
galaxies (e.g., Zhang etal. 2012). Use of the FB for iden- 


tifying X-ray point sources is appropriate for our early-type 
galaxy sample, since these galaxies have negligible contam- 
ination from hot gas (see, e.g.. Fig. 1) and little absorption 
at ^ 1 keV from cold interstellar gas and dust. To com- 
pute XLFs for each of our galaxies, we made use of the FB- 
detected sources in the main Chandra catalogs presented in 
§3.1. Spatial variations in X-ray sensitivity, due to variations 
in local backgrounds, PSFs, effective exposures (e.g., chip 
gaps and bad pixels and columns), and source crowding, com- 
bined with incompleteness can have significant effects on the 
shape of the XLFs at X-ray luminosities that are within a fac- 
tor of ~ 10 of the detection limits. These effects need to be 
accounted for to properly compute the XLFs. Our approach, 
described in detail below, is to fit the observed XLFs with 
models of the intrinsic LMXB XLFs convolved with com- 
pleteness functions derived from simulations. 

5.1. X-ray Background and Sensitivity Maps 

Our first step in computing completeness functions was to 
construct FB background and sensitivity maps for each of the 
three galaxies in our sample. Background maps were con- 
structed by first masking out all X-ray point sources within the 
merged FB images using circular masking regions of size 1.1 
times the radius that encloses 99% of the PSF. These mask- 
ing regions were computed and utilized by ae in the process of 
computing local background estimates for the point sources. 
After masking out the point sources, we filled the masked re- 
gions with Poisson noise that was characteristic of the local 
backgrounds as computed by ae. The background maps en- 
abled us to create sensitivity maps, which contain pixel-by- 
pixel estimates of the minimum number of counts ,v mm that a 
hypothetical source would need to be detected given the lo- 
cal background. To compute ,v mm at each location in a given 
image, it is necessary to estimate local values of /? src and b ext 
(via Eqn. (1)) that are similar to those that would be estimated 
following the ae approach used for the X-ray detected sources. 
Unfortunately, the ae procedure for extracting these values at 
every location in the image is computationally prohibitive, 
since source photometry is performed on an observation-by- 
observation basis using computationally intensive ray-tracing 
algorithms to generate local PSFs and approximations to the 
90% EEF PSF contours. To overcome this issue, we estimated 
/Arc by extracting background counts from the background 
maps using circular apertures with sizes that encompass the 
90% EEF for a point source. Appropriate values of b ext were 
estimated at each location (pixel) by extracting the maximum 
value of b ext from the main Chandra catalogs for sources with 
off-axis angles 6 = 8 P ± 0(25, where 9 p is the pixel off-axis an- 
gle. Given values of Z> src and /? ext , we numerically solved the 
relation P(x > s m i n ) = 0.004 (i.e., Eqn. (1)) to obtain .s m ; n at 
each image location. In this manner, we constructed a spatial 
sensitivity map consisting of s m j n values. 

5.2. X-ray Completeness Functions 

The above maps allow us to determine spatially whether a 
source would have been detected for a given number of total 
counts. However, these maps alone cannot directly be used to 
assess source completeness and encapsulate subtleties related 
to our complex selection process (i.e., wavdetect plus ae 
pruning) and source confusion in crowded regions. To as- 
sess better source detection completeness as a function of net 
source counts for each Chandra data set, we made use of sim- 
ulated data sets and a Monte Carlo procedure that replicated 
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Figure 4. Fraction of recovered simulated sources as a function of net source counts for our three early-type galaxies (filled circles with 1 a error bars). A 
simulated source was considered “recovered” if (1) the total number of on-source counts satisfied the binomial probability selection criterion of equation 1, (2) 
the source was detected by wavdetect using a false-positive probablity threshold of 10~ 5 , and (3) the source had net counts that were within 3cr of the input 
counts (required to identify confused sources). Each color represents a different offset range from the nuclear position, as annotated in each plot. Our best-fit 
models have been plotted as dashed curves. For a complete description of our simulations, see § 4.2. 

non-negligible effect on the source confusion properties of the 
original images and allows for 35 times more efficient execu- 
tion of the simulation process. For each X-ray source added 
to a given mock image, we assigned random right ascension 
and declination values that were probabilistically offset from 
the galaxy nucleus in accordance with the offset distribution 
of the main Chandra catalog. Each X-ray source was given 
a random number of counts between Ril-100 counts. The 
counts from each source were added to the simulated image 
canvases using the mar>Fl ( version 4.5) ray-tracing code. In 
this procedure, the PSF of each source was modeled assum- 
ing a roll angle and aim point that were selected from one of 
the galaxy ObsIDs (Table 2). Selection of the roll angle and 
aim point for each source was done probabilistically with the 
probability of selection being directly proportional to the ex- 
posure time of each of the observations. 

For each galaxy, we performed circular aperture photome- 
try on the 28,000 mock X-ray sources (i.e., 35 sources in each 
of the 800 mock images) using a circular aperture with a ra- 
dius encompassing «90% of the encircled energy. The num- 
ber of counts measured for each source was then compared 
with the number of counts needed for a source detection from 
the sensitivity maps to see if it satisfied the detection criterion 
adopted in equation (1). We then searched each of the 800 
mock images using wavdetect at a false-positive probabil- 
ity threshold of 10~ 5 and made wavdetect source catalogs 
for each image. We considered a simulated source to be “re- 
covered” if (1) it was detected in the wavdetect runs, (2) 
the total number of counts exceeded the value required for a 
detection based on the sensitivity maps, and (3) the net counts 
derived were within 3 a of the known input counts. Crite- 
rion 3 was adopted to account for source confusion, since 
sources that are significantly confused with brighter sources 
in the original image will have an excess of source counts due 
to the unrelated source. We consider such sources to not be 
recovered. 



S = 0.5-8 keV Flux (ergs cm 2 s ') 


Figure 5. 0.5-8 keV number counts measurement for X-ray sources within 
the galactic footprints of our early-type galaxies that were optically classi- 
fied as unrelated background objects ( see annotations ). For comparison, we 
show the expected number counts of unrelated background objects ( dotted 
curve), and lcr error envelope ( grey region ), based on number counts from 
extragalactic X-ray surveys (Georgakakis et al. 2008). This comparison indi- 
cates that we are able to identify the majority of unrelated background X-ray 
sources using our methods. The deficit of number counts in NGC 3379 may 
be due to cosmic variance (see discussion in §5.3). 


the source detection and pruning processes that were used to 
create our source catalogs. 

We began by generating 800 mock FB images for each 
galaxy with 35 X-ray sources artificially added to the origi- 
nal FB images. Ideally, we would have generated more mock 
images with only 1 X-ray source implanted in each mock im- 
age, so that the confusion properties would not be changed by 
the addition of sources. However, we found empirically that 
the addition of 35 sources per image both avoids having any 


1 See http://space.mit.edu/cxc/marx/ for marx simulator details 
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Figure 6. Cumulative observed X-ray luminosity functions (N(> Lx); thick black curves ) for the GC plus field LMXB populations of the three early-type 
galaxies in our sample. Each observed XLF (gray circles with la Poisson error bars) includes X-ray detected point sources within the extent of the galaxy after 
excluding unrelated background sources and sources within the central nuclear regions (see text). The XLFs were fit using an intrinsic XLF model convolved 
with our completeness functions following the procedure described in §5.4. 


Our simulations revealed that the central «10"-15" of the 
three early-type galaxies suffer catastrophic source confusion 
due to a high density of bright X-ray sources in the nuclear 
regions. As discussed in §3.2, these regions were consid- 
ered problematic from an optical perspective, since the bright 
emission associated with the stellar fields made it impossible 
to identify potential GC and unrelated background candidates. 
These regions are problematic from the X-ray perspective be- 
cause (1) it is not possible to recover accurately the true X-ray 
luminosity distributions of the sources in these regions; (2) 
the X-ray point sources may be formed via dynamical inter- 
actions in the high stellar density environments (e.g., Voss & 
Gilfanov 2007; Zhang etal. 2011), which are a different for- 
mation channel from the field LMXBs that we are interested 
in studying here; and (3) among the X-ray point sources of a 
given galaxy may be a low-luminosity AGN that could con- 
taminate the LMXB XLFs at the high-luminosity ends. These 
considerations further support our choice to exclude these re- 
gions from our XLF computations below. 

In Figure 4, we show the recovery fraction as a function of 
net source counts in radial bins for our three early-type galax- 
ies. This plot was generated using circular annuli including 
only sources that were within the elliptical regions defined in 
Table 1; for the reasons discussed above, the exclusion re- 
gion defined in Column (8) of Table 1 has not been shown. 
For NGC 3115 and 3379, we find that source recovery be- 
comes more complete moving from the innermost annulus to 
the second innermost annulus. This is due to confusion be- 
coming less severe at larger radii. The completeness again 
deteriorates going to even larger radii as the size of the PSF 
increases somewhat with radial distance (N.B., the exposure- 
weighted Chandra aimpoints are all within «5" of the galaxy 
centers). For NGC 3384, confusion is not a strong issue for 
the innermost annulus; therefore, the completeness simply de- 
clines with radius. 

5.3. Unrelated Background X-ray Sources 

One of the key strategies in our XLF computations is 
to identify directly and eliminate from consideration X-ray 


point sources that are associated with unrelated background 
objects. Direct optical identification of background X-ray 
sources is preferred over statistical methods, in which a back- 
ground component is simply subtracted from the total XLFs, 
since the galaxies subtend only small solid angles on the sky 
( ^ 0.01 deg 2 ) and are subject to large variations in their back- 
ground populations due to cosmic variance (e.g., Moretti et al. 
2009). To test the efficacy of our background identifica- 
tions, we constructed the X-ray number counts (log/V-log.S - ; 
the number of X-ray sources per unit area N with X-ray flux 
brighter than S) relations for each of the three galaxies (and 
the combined sample) using the X-ray sources classified as 
unrelated background objects (see §4). For a given galaxy, 
we limited this calculation to a flux regime in which we could 
have detected sources across «75% of the galactic extent. 

Figure 5 shows the log/V-log.S’ for the unrelated back- 
ground sources in the early-type galaxy footprints (as well 
as the total galaxy sample) compared with expectation based 
on extragalactic blank-field X-ray surveys (e.g., Georgakakis 
etal. 2008). It is apparent that the number counts for un- 
related background sources for NGC 3115 and 3384 are in 
good agreement with the expectation, while the NGC 3379 
number counts fall a factor of «2.5-3 times below the 
expectation across the full X-ray flux range. We ex- 
pect that our optical identification and classification pro- 
cess will be most complete and accurate at bright fluxes 
(/0.5-8 kev ~ few x 10~ 15 erg cm -2 s -1 ) where unrelated back- 
ground counterparts are likely to be relatively unobscured 
AGN with bright optical counterparts (e.g., Steffen etal. 
2006; Xue etal. 2011). Since the deficit of X-ray sources 
in NGC 3379 is apparent at these high fluxes, it is likely that 
the deficit is real. As noted in §3.3, NGC 3379 also has a 
factor of «2 deficit in the source density of detected optical 
sources, lending strong support to this conclusion. Given the 
areal coverage of NGC 3379, we would expect cosmic vari- 
ance at a level of ~ 3 0-3 5% (see Moretti etal. 2009). There- 
fore the factor of >, 2 deficit for NGC 3379, if due to cosmic 
variance, would represent a 3<r deviation from the expected 
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value (assuming the area-dependent cosmic variance distribu- 
tion is symmetric and Gaussian). 

The above number counts analysis suggests that the ma- 
jority of unrelated X-ray background sources are likely to be 
classified as such by our procedure. In our computations of 
LMXB XLFs, described below, we remove unrelated back- 
ground X-ray sources and do not try to correct for remaining 
contributions from unrecovered background X-ray sources. In 
the faint X-ray flux regime, where our recovery process is ex- 
pected to be the least complete, background X-ray sources not 
classified as such are greatly outnumbered by field LMXBs 
and therefore unlikely to have material effects on our results. 


5.4. X-ray Luminosity Function Fitting 

For each galaxy, we constructed observed XLFs using the 
FB-detected sources in the main Chandra catalogs that were 
within the galactic footprints (excluding sources in the near 
vicinities of the galactic nuclei; see Column (8) of Table 1). 
For each galaxy, we chose to construct and study the XFFs for 
three FMXB population categories, including GC plus field 
(total) FMXBs, field LMXBs, and GC LMXBs. However, 
for NGC 3384, we found only three GC LMXB candidates, 
which were all undetected in the FB, and therefore did not 
attempt to construct their GC LMXB XLF. In all cases, we 
excluded sources with optical counterparts classified as unre- 
lated background objects. For illustration, we have displayed 
the cumulative observed XLFs (i.e., N(> L x )) for the GC plus 
LMXB populations in Figure 6 (gray points with error bars). 

For each of the three LMXB sub-populations, we modeled 
intrinsic XLFs using a broken power-law form that followed 


dN / dLx = 


K(L x /L lef T a ' (L x < L b ) 

X(Lb/L re f)“ 2_Ql (L x /Lrcf) _Q!2 , (L x > Lb) 


where K is the normalization at reference luminosity L re f = 
10 38 erg s _1 , «i is the faint-end slope. Lb is the break lumi- 
nosity, and a 2 is the bright-end slope. 

We chose to model the observed XLFs using a forward- 
fitting approach, in which we convolved an appropriately 
weighted completeness function with our model XLFs before 
fitting to the observed XLFs. As demonstrated in § 5.2 and 
Figure 4, our X-ray completeness to sources depends upon 
their offsets from the galactic centers. For a given galaxy, we 
therefore statistically weighted the contributions from model 
XLFs at each annulus according to the observed distributions 
of X-ray point sources to obtain the galaxy-wide weighted 
completeness function £(L X ). Here we utilized only luminous 
X-ray sources in regimes where our point source populations 
are highly complete at all offsets from the galactic centers to 
determine the annular weightings. Formally, we computed 
£(L X ) using the following relation: 


£(L X ) = ^/ re 


v,i(Fx) x W„ 


( 2 ) 


where / reC ov,iCL x ) is the recovery-fraction curve for the ;th an- 
nular bin and w, is the fraction of total number of galaxy-wide 
sources within the ;th annuluar bin based on the observed 
point-source distributions. 

Using the sherpa modeling and fitting package available 
through ciao, we modeled each observed XLF using a mul- 
tiplicative model 

dN/dL x ( observed) = £(L x )dN /dL x (model). (3) 


For a given subsample of sources, we constructed the ob- 
served dN/dLx( observed) using a very small constant bin 
of SL x spanning the minimum luminosity of the subsample 
to 10 times the maxmimum luminosity. Therefore the vast 
majority of bins contained zero sources up to a maximum 
of three sources. We evaluated the goodness of fit for our 
double power-law models using the Cash statistic (cstat; 
Cash 1979). We tested to see if in any cases a single power- 
law provided a statistically superior fit over the broken power- 
law, but were unable to find any such cases. We did find that 
for the field LMXB XLF in NGC 3115, the best-fitting broken 
power-law produced an effectively infinite high-luminosity 
slope ( 0 : 2 ) beyond the most luminous detected source. For 
this subclass, we utilized a single power-law with a high- 
luminosity cut-off (see Table 5). In Figure 6, we display the 
observed XLFs and best-fitting observed and intrinsic XLF 
models for the GC plus field LMXB populations for our three 
early-type galaxies. In Table 5, we tabulate the best-fit param- 
eters and 1 a errors for fits to the GC plus field LMXBs, field 
LMXBs, and GC LMXBs. 

6. RESULTS 

In Figure 7, we show the L x -normalized cumulative LMXB 
XLFs (N(> L x )/Lk) for our early-type galaxies for GC plus 
field and field LMXB populations. The XLFs displayed in 
Figure 7 have been corrected for completeness, near the lu- 
minosity limits of each galaxy, using the completeness cor- 
rections derived from our best-fitting XLFs (see §5.4). We 
find that the GC plus field LMXB XLFs (Fig. la) for all 
three galaxies appear to be similar in normalization and shape 
across L x w 10 37 — 10 39 erg s -1 . These estimates are in good 
agreement with the best-fit XLFs derived by Zhang etal. 
(2012) for a sample of 20 massive early-type galaxies (gray 
curve). However, when we remove the contributions from GC 
LMXBs, we observe significant galaxy-to-galaxy variations 
in the resulting field LMXB XLFs. 

In Figure lb. we display the field LMXB XLF. We find 
that the shapes and normalizations of the field LMXB XLFs 
are similar for the two old early-type galaxies, NGC 3115 
and 3379, with some differences becoming apparent at 
L x < (2^1) xlO 37 erg s _l . For both of the old early-type 
galaxies, we observe sharp declines in the field LMXB XLFs 
above L x ~ (0.5-1) xlO 38 erg s -1 . By contrast, the field 
LMXB XLF for the young early-type galaxy NGC 3384 is 
generally elevated compared to those of the old early-type 
galaxies and includes an extension to higher X-ray luminosi- 
ties (Lx ~ 7 x 10 38 erg s _1 ); however, as is apparent from Fig- 
ure lb, the cumulative field XLF of NGC 3384 approaches 
that of NGC 3379 at L x ss 3 x 10 37 erg s _I . A similar excess 
of luminous (L x > 5 x 10 38 erg s -1 ) LMXBs in young early- 
type galaxies was previously observed by Kim & Fabbiano 
(2010) and Zhang etal. (2012) using GC plus field LMXB 
populations. Our analysis shows that this excess may appar- 
ently be attributed to a strong excess in field LMXBs, and the 
removal of GC LMXBs seems to amplify this excess of bright 
X-ray sources. This result clearly illustrates that HST data are 
crucially needed to identify directly LMXB counterparts and 
provide meaningful tests of population synthesi models. By 
selection, the galaxies in our sample have relatively low Sn 
(Sn ~ 3), so the large elliptical galaxies that have been stud- 
ied well thus far will contain even larger contributions from 
GC LMXBs. 

In Figure 8, we display the GC LMXB cumulative XLFs 



Dependence of Field LMXB XLFs with Stellar Age 


13 


Table 5 

Best-Fit Parameters to X-ray Luminosity Functions (XLFs) 


GC + Field (Total) LMXBs Field LMXBs GC LMXBs 


Galaxy K cx.\ a.2 L K ot\ oii L b K ol\ 0L2 L b 

NGC3115 6.8±1.B 1.5±0.1 2.2d=0.7 1.76±0.02 4.6±1.2 1.5±0.1 ~ 1.853^;^* 2.8±1.2 1.3±0.2 1.7±0.6 1.76=1=0.19 

NGC3379 9.4±3.2 1.6±0.1 2.3±0.4 0.60=1=0.19 7.1±2.1 1. 7=1=0. 1 3.3=fcl.l 0.84=1=0.17 4.6±5.1 0.6±0.4 1.7±0.4 0.41=1=0.26 

NGC3384 6.1±3.1 1.3±0.4 2.0±0.4 0.68±0.23 6.1±3.1 1.3±0.4 2.0±0.4 0.683=0.23 


NOTE — K and Lj are in units of (10 38 erg s 1 ) 1 and 10 38 erg s 1 . respectively. All errors are 1 a and are based on the fitting procedure described in §5.4. For 
NGC 3384, all FB detected LMXBs have been classified as field candidates; therefore the GC plus field and field only fit parameters are identical. 

*The best-fit model to the field LMXB XLF for NGC 3115 was a single power-law with normalization K , slope ol\, and high-luminosity cut at the quoted Lb 
value. 
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Figure 7. A'-band-luminosity-normalized cumulative XLFs, N(> Lx)/Lk, for our early-type galaxies in two LMXB population categories: (a) GC plus field 
LMXBs and ( b ) field LMXBs. Each XLF has been summed down to an X-ray limit that is roughly 50% complete, and the displayed XLFs have been corrected 
for completeness following the methods described in §5.4. The three early-type galaxies have GC plus field (total) LMXB XLFs consistent with each other and 
previous results (Zhang etal. 2012). However, when looking at only the field LMXB XLF contribution, we find that the young early-type galaxy NGC 3384 has 
a more numerous and luminous field LMXB population than the old early-type galaxies NGC 3115 and 3379 (see §6 for more details). 
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Figure 8. GC LMXB cumulative XLFs normalized by Lk and Sn ( col- 
ored curves) with field LMXB cumulative XLFs, normalized by Lk only, 
displayed for comparison {gray curves). Linestyles and colors follow the 
same scheme as in Figure 7, with the exception of NGC 3384, for which we 
only obtain an upper limit on the GC cumulative XLF at the detection limit 
(filled blue circle with arrow). We find that when normalized by Lk and Sn, 
the GC LMXB XLFs are all consistent with eachother in shape and normal- 
ization and have overall slopes shallower than the field LMXB XLFs. 


normalized by Lk and Sn (N(> Lx)/Lk/Sn)- We note that 
NGC 3384 had no GC LMXBs detected in the FB; we there- 
fore show only the 3er upper limit to its GC XLF at the de- 


tection limit. From Figure 8, we see that the GC LMXB 
XLFs for NGC 3115 and 3379 appear to be broadly consis- 
tent in their shape and normalization; the upper limit on the 
GC LMXB XLF for NGC 3384 is consistent with the other 
galaxies. For the purpose of comparing XLF shapes, we have 
plotted in light gray N(> Lx)/Lk (without the Sn normaliza- 
tion) for the field LMXBs. We find that the GC LMXB XLFs 
of NGC 3115 and 3379 have significantly shallower slopes 
and extend to higher X-ray luminosities than the field LMXB 
XLFs. This result is consistent with previous results that com- 
pared GC and field LMXB XLFs (e.g., Kim et al. 2009; Voss 
etal. 2009; Paolillo etal. 2011). 

The above findings indicate that the young early-type 
galaxy NGC 3384 contains a more numerous and X-ray lumi- 
nous population of field LMXBs per unit K-band luminosity 
(stellar mass) than the older early-type galaxies NGC 3115 
and 3379 (see §6). Although the integrated X-ray luminos- 
ity per unit K-band luminosity Lx/Lk can be measured to 
high precision for the LMXB populations in the three galax- 
ies in our sample, and our measurements clearly indicate that 
the field LMXB (Lx/Lk ) geld is larger for NGC 3384 than for 
NGC 3115 and 3379, it is important to account for the fact that 
the Lx value for a given galaxy can be substantially influenced 
by the brightest few sources. If we consider the distribution 
of point-source luminosities as a realization of an underlying 
XLF, then we might expect large galaxy-to-galaxy variations 
in Lx/Lk simply due to “statistical noise” (see Gilfanov et al. 
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Figure 9. ( a ) Integrated LMXB X-ray luminosity per Lk versus stellar age for the field LMXBs, (Lx/L/Ofieid* in the three galaxies in our sample {filled circles 
with annotations ) and the total (GC plus field) LMXBs, (Lx/LkX otab in the Boroson et al. (2011) sample ( gray squares. The three filled gray squares indicate the 
Boroson etal. (2011) (Lx/LkXo tai values for the three galaxies in our sample. The dashed black curve and surrounding gray region is the expected relation and 
spread from the population synthesis models of Fragos etal. (2013b), and the solid curve shows the best-fit version of this curve renormalized to our data. The 
dotted gray curve is the best-fit relation to the GC plus field LMXB (Lx/L/cXotai versus age data (i.e., the gray squares) and the long-dashed curve shows how this 
correlation would change after subtracting the expected contributions from GC LMXBs (see Eqn. (7) and discussion in § 7 for details). ( b ) The observed relation 
between GC Sn and stellar age with the best-fitting model {dotted gray curve). We find that the apparent correlation in Figure 9 a between {Lx/Lk) total and stellar 
age can be fully explained by contributions from GC LMXBs and the underlying correlation between GC Sn and stellar age shown in Figure 9b. Furthemore, the 
contributions from GC LMXBs could mask an underlying anticorrelation between {Lx/Lk ) field and stellar age (e.g., the long-dashed curve in Figure 9a). 


Table 6 

Integrated X-ray Luminosities 


log Lx log Lx/Lk 

Galaxy (log ergs s -1 ) (log ergs s -1 L^U) 

GC + Field (Total) 

NGC 31 15 39.58+/- 0.19 28.90+/- 0.19 

NGC 3379 39.51 +/- 0.12 28.75+/- 0.12 

NGC 3384 39.37 +/- 0.26 28.87 +/- 0.26 

Field 

NGC 3115 39.08 +/— 0.10 28.40+/-0.10 

NGC 3379 39.35 +/- 0.07 28.58 +/- 0.07 

NGC 3384 39.37 +/- 0.26 28.87 +/- 0.26 

GC 

NGC 3115 39.52 +/- 0.34 28.84 +/- 0.34 

NGC 3379 39. 17 +/- 0.50 28.40 +/- 0.50 


2004 for a detailed discussion on this issue). For example, a 
large sample of galaxies with identical underlying XLFs will 
contain a spread of Lx values due to this statistical noise. 

To quantify how this statistical noise influences our results, 
we made use of a Monte Carlo technique to measure the ex- 
pected standard deviation of Lx/Lk for a given XLF. For each 
of our LMXB populations in each galaxy, we constructed 
20,000 simulated XLFs that contained distributions consis- 
tent with the best-fitting XLF model (see §5.4 and Table 5). 
For a given simulation, we first computed the expected num- 
ber of LMXBs that would have Lx 10 36 erg s -1 and per- 
turbed this value using Poisson statistics. We then assigned 
Lx values to each of the sources probabilistically using the 
best-fitting XLF model as our underlying probability distribu- 
tion function. Finally, the sources within a given simulation 
were summed to give an integrated Lx value. The resulting 
distribution of 20,000 Lx values provided us with the effects 
of the statistical noise on our results. In Table 6, we have tab- 
ulated the integrated Lx and Lx/Lk values and their errors. 


which include the statistical noise, and in Figure 9a, we have 
plotted the field LMXB (Lx/Lk) field versus stellar age (and ler 
errors) for the galaxies in our sample. We have found that the 
resulting errors on log Lx values are well described by Gaus- 
sian functions, implying asymmetric Lx errors. Considering 
these errors, and following the error propogation methods de- 
scribed §2.4.5 of Eadie etal. (1971), we find that the field 
LMXB (Lx/Lk ) field for NGC 3384 is a factor of 2.0 and 
3.0 /j 4 times higher than NGC 3115 and 3379, respectively. 

7. DISCUSSION 

Our results support a scenario in which the field LMXB 
X-ray luminosity per unit X-band luminosity (Lx/L k )i\m de- 
clines with increasing stellar age, as predicted by population 
synthesis models (filled circles in Fig. 9a). In Figure 9 a we 
have displayed the expected mass-weighted stellar age depen- 
dent (Lx/Lk) field from the population synthesis models from 
Fragos etal. (2013b) (dashed curve with gray shaded region). 
This curve was derived using equation (4) of Fragos etal. 
(2013b) and an assumed M/Lk « 0.66 Mq/Lk,q, consistent 
with our galaxies. We find that the population synthesis mod- 
els predict that galaxies with mass-weighted stellar ages of 
«3.2 Gyr will have (Lx/Lk) field values that are a factor of «3 
times higher than «8.5 Gyr galaxies, which is consistent with 
our finding. We note that the absolute scaling of our field 
LMXB (Lx/Lx)fieid values appear to be a factor of ~ 1 .5- 
2 times lower than the maximum likelihood model adopted 
in the Fragos models. However, it is important to note that 
these models were calibrated to match local total (i.e., GC 
plus field) (Lx/Lk) total LMXB scaling relations, which as we 
have shown in this study, are affected somewhat by the GC 
LMXB populations. Furthermore, the range of statistically 
acceptable population synthesis models (i.e., the shaded re- 
gion in Fig. 9a; see Fragos et al. 2013b for details) covers our 
field LMXB (Lx/Lk) field values. Using the three (Lx/L^)fl e id 
values derived from this study, we adjusted the normalization 
term in equation (4) of Fragos et al. (2013b) to fit our results. 
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Using x 2 minimization, we found that a factor of ss 1.6 down- 
ward shift of the normalization would provide a reasonable fit 
to our data ( solid curve in Fig. 9a). 

Our pilot study provides suggestive evidence that the X-ray 
emission from field LMXB populations declines with stellar 
age; however, it is absolutely critical to expand this study to a 
larger sample of early-type galaxies using comparable Chan- 
dra and HST imaging before we draw statistically robust con- 
clusions. As we noted in § 1, previous studies have found that 
the total (GC plus field) (Lx/Lk) total increases with stellar age. 
In Figure 9a, we have displayed the (Lx/Lk ) total versus age for 
the Boroson et al. (201 1) sample. It is clear that indeed there 
is a positive correlation between (Lx/Lk \ ota i and stellar age; 
a Spearman rank test indicates that the correlation is signifi- 
cant at the > 99.9% confidence level. However, as pointed out 
by Zhang etal. (2012), this correlation is strongly influenced 
by a secondary correlation of GC Sn with stellar age (see Fig- 
ure 9b). As discussed in §6, we have shown that the Lk and Sn 
normalized GC XLFs for the galaxies in our sample are con- 
sistent with each other, suggesting that the GC (Lx/Lk/Sn ) gc 
may be a “universal” quantity for early-type galaxies gen- 
erally. A similar conclusion was reached by Irwin (2005), 
who determined the fi-band luminosity normalized quantity 
(Zfl. 3 -io kev/Cg/ Sn)gc = (2.72 ±0.47) x 10 29 erg s _1 L]/ Q for a 
sample of 12 galaxies. If we assume Lx/Lo. 3-10 kev ~ 0.8 and 
Lb/Lk w 0.18 Lb,q/Lk,q, appropriate for our X-ray binaries 
and galaxies, respectively, we can convert the Irwin (2005) re- 
lation to (Lx/L k /Sn)gc = (3.70 ±0.64) x 10 28 erg s -1 Lj^g. 
Using the data provided in Tables 1 and 6, we find that 
(Lx/L k /S n ) gc « 3.0^ x 10 28 erg s -1 L~ K \ Q and 1.9^; 2 x 
10 28 erg s -1 LTf ) Q for NGC 3115 and 3379, respectively, both 
formally consistent with the Irwin (2005) estimate. 

Given the measurements of GC LMXB (Lx/Lk/Sn)gc and 
correlations between GC plus field (Lx/Lk ) total and .S'v with 
stellar age, we can roughly estimate how the field (Lx/Lx)fi e id 
would be expected to scale with stellar age. This can be done 
by solving the following system of equations: 


1 Og( /.\ L/, itotal ~ A + B Lge? 

(4) 

(Lx/Lk) gc ~ C Sn, 

(5) 

Sn r '^ L) t a g e , 

(6) 

(Lx/Lk) field = (Lx/Lk\o\a-(Lx/Lk)gCi 


-A (Lx/Lk) fie id ~ 10 4+/f ** -CxD f age , 

(7) 


where A, B, C, and D are constants. Using the Boroson etal. 
(2011) data presented in Figures 9 a and 9b, we determined 
A ss 28.84 log erg s _1 Lj^ 0 , B « 0.032 log erg s _l L]/ & Gyr -1 , 

and D ss 0.42 Gyr -1 . The dotted gray curves in Figures 9 a 
and 9b illustrate the best-fit solutions to equations (4) and (6), 
respectively. If we adopt C « 3 x 10 28 erg s _1 L^ Q (see above 
for justification) and insert this set of values into equation (7), 
then we obtain an estimate for (L X /L K ) field- which we show 
as the long-dashed curve in Figure 9a. This exercise illus- 
trates that the positive correlation of (Lx/Lk ) total with stellar 
age is dominated by contributions from GC LMXBs and that 
it is possible to mask an underlying negative correlation be- 
tween (Lx / Lk) H eld and stellar age. We note, however, that 
the above calculation is very sensitive to the set of constants 
that are adopted. Although each of the constants above were 
derived directly from the data, it is possible to perturb a con- 
stant in a manner entirely consistent with the data and obtain 
a qualitatively different result. For example, if we adopted 


C = 2 x 10 28 erg s _1 L^' 0 , the value obtained for NGC 3379, 
we would find the long-dashed curve in Figure 9 a to have 
a flat slope. Or if we adopted C = 4 x 10 28 erg s' 1 Lk\q, 
a value higher than but consistent with those obtained for 
NGC 3115 and 3379 and the Irwin (2005) result, the re- 
sulting (Lx/ Lx) field from equation (7) would be negative at 
f a ge > 6.7 Gyr. Therefore, the above calculation should be 
taken as a proof of concept and not a well constrained mea- 
surement of the stellar age dependence of (Lx/Lk) field- To 
obtain an accurate measurement of how (Lx/Lk) field varies 
with stellar age will require direct measurements of the field 
LMXB populations in a larger sample of galaxies using the 
combination of HST and Chandra data sets. Obtaining such 
an accurate measurement will be important for constraining 
X-ray binary population synthesis models. For example, the 
above inferences for how (Lx/Lk) field changes with stellar age 
(i.e., the long-dashed curve in Fig. 9) appears flatter than the 
current X-ray binary population synthesis predictions at the 
youngest stellar ages. Direct measurements of (L\/ Lk)um 
for several galaxies are required to show whether this is in- 
deed the case and whether these predictions will need to be 
revised. 

8. SUMMARY 

We have utilized deep Chandra and HST observations of 
the three nearby early-type galaxies NGC 3115, 3379, and 
3384 to identify and study subpopulations of LMXBs that are 
found in the galactic field environments. Our galaxy sample 
was selected to span a large range of stellar age («3-10 Gyr) 
so that we could measure how the emission from field LMXBs 
depends on stellar population ages. Our results can be sum- 
marized as follows; 

• We find that LMXBs formed dynamically in GCs and 
unrelated background X-ray sources can significantly 
impact the overall XLF shapes and X-ray power out- 
put from early-type galaxies, even in low-5^ galaxies 
like the three sources studied here (Sn ^ 3). This find- 
ing indicates that for studies of field LMXBs, it is crit- 
ical to remove GC LMXBs and unrelated background 
sources using the combination of HST and Chandra 
data. This is especially the case for the majority of 
early-type galaxies that have been studied in the liter- 
ature, which have Sn ~ 2 (see Fig. 7 and §§4 and 6). 

• After normalizing for stellar mass (via the L'-band lumi- 
nosity), we find that the field LMXB population in the 
young («3.2 Gyr) galaxy NGC 3384 is more numer- 
ous and extends to higher luminosities than the older 
(«8.5 Gyr) galaxies NGC 3115 and 3379. This re- 
sult is consistent with previous studies that have found 
a similar excess of luminous X-ray sources in young 
( ^ 6 Gyr) versus old ( X, 6 Gyr) early-type galaxies 
(e.g., Kim & Fabbiano 2010; Zhang etal. 2012) (see 
Fig. lb, and §6). 

• Using the field LMXB populations in our galaxies, we 
measure the integrated (Lx/Lx)fi e id for NGC 3384 to 
be a factor of «2-3 times higher than that found in 
NGC 3115 and 3379. This result is consistent with 
the trends expected from the X-ray binary population 
synthesis models published by Fragos etal. (2013a, 
2013b), which predict that the (Lx/Lx)fi e id will be a fac- 
tor of «3 times higher for galaxies with stellar ages of 
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~3 Gyr compared to those with ss9 Gyr (see Fig. 9 a 
and §§6 and 7). 

• We find that the GC LMXB XLFs for NGC 3115 and 
3379 are consistent with each other when normalized 
by L k and .S' v . We did not detect enough GC LMXBs 
in NGC 3384 to study the shape of its GC LMXB 
XLF; given its low Sn, our lack of GC LMXB detec- 
tions is not surprising. We constrained the mean GC 
LMXB (L x /L k /S n )gc « (2-3) x 10 28 erg s" 1 L K \ Q for 
these galaxies, in good agreement with past studies (Ir- 
win 2005) (see Fig. 8 and §§6 and 7). 

• Previous studies have found that the total GC plus 
field LMXB (L X /L K ) tota i was positively correlated with 
stellar age, which was considered an apparent con- 
flict with population synthesis predictions. This ob- 
served correlation is largely driven by an underlying 
correlation between Sn and stellar age, which masks 
any underlying correlations of (L x /L x ) field with stel- 
lar age. We showed that if we consider our con- 
straint on (L x /L x /Sn)gc and the .S'y-ste 1 1 ar age corre- 
lation, it is plausible that (L x /L K ), lM = (L X /L K ) tot ai - 
(L x /Lk/Sn)gc declines with increasing stellar age. 
However, a robust measurement of (L x /L x ) field versus 
age will only come from a detailed Chandra and HST 
study of the field LMXB populations from a larger sam- 
ple of early-type galaxies (see Fig. 9 and §§1,6, and 7). 
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